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ABSTRACT 

We study possible anomalous top-quark couplings generated by SU(2) x U(l) 
gauge-invariant dimension-6 effective operators, using the final 6-quark momentum 
distribution in 77 —>■ ti — > bX. Taking into account non-standard tij, tbW and 
77-ff couplings, we perform an optimal-observable analysis in order to estimate the 
precision for the determination of all relevant non-standard couplings. 
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1. Introduction 

Linear colliders of e + e~ are expected to work as top-quark factories, and therefore a 
lot of attention has been paid to study possible non-standard top-quark interactions 
through ee — > ti (see, for instance, P] I2| and their reference lists). An interesting 
option for such e + e~ machines could be that of photon-photon collisions, where 
initial energetic photons are produced through electron and laser-light backward 
scatterings jSHU- 

This type of colliders presents remarkable advantages for the study of CP vi- 
olation. In the case of ee collisions, the only initial states that are relevant are 
CP-even states ^l/rGr/l) under the usual assumption that the electron mass 
can be neglected and that the leading contributions to ti production come from 
s-channel vector-boson exchanges. Therefore, all CP-violating observables must 
be constructed from final-particle momenta/polarizations. In contrast, a 77 col- 
lider offers a unique possibility of preparing the polarization of the incident-photon 
beams, which can be used to construct CP-violating asymmetries without relying 
on final-state information. 

This is why a number of authors have considered top-quark production and 
decays in 77 collisions in order to study i) Higgs-boson couplings to the top quark 
and photon j^-JH], or ii) anomalous top-quark couplings to the photon [T2]-|14j. 
However, what is supposed to be observed in real experiments is combined signals 
that originate both from the process of top-quark production and, in addition, 
from its decays. Therefore, in our latest article [TS] we considered 77 — > ti — > £ + X, 
including all possible non-standard interactions together (production and decay), 
and performed a comprehensive analysis as model-independently as possible within 
the effective-Lagrangian framework of Buchmiiller and Wyler 16 . 

In this letter, we will carry out an optimal-observable (00) analysis, using the 
final 6-quark momentum distribution, as a complementary work to What we 
have to do for this purpose is similar to what has been done in However, in 
the case of the bX final state, we can expect to obtain independent and valuable 
information since there is no branching-ratio suppression for t —>■ bW, in contrast 
to the analysis with the final lepton. One might say that using the final 6-quark 



distribution is not that effective, since the determination of the 6-quark momentum 
is more challenging than that of charged leptons. However, in any case, it is crucial 
to tag the final 6-quark efficiently in order to distinguish the top-quark production 
from the main background of W + W~ production That is, we cannot study 
top-quark events without good information on the final 6-quark, which makes our 
analysis realistic. 

2. Framework 

We use the effective low-energy Lagrangian [TBI EE] to describe possible new-physics 
effects. Following this approach, we consider the Standard- Model (SM) Lagrangian 
modified by the addition of a series of SU(2) x (7(1) gauge-invariant operators (9, 
whose coefficients parameterize the low-energy effects of the underlying high-scale 
physics. 

Since the detailed description of this framework was presented in we only 
mention here that the largest contribution comes from dimension-6 operators, and 
that these lead to the following Feynman rules for on-shell photons, which are 
necessary for our calculations: 

(1) CP-conserving ti r ) vertex 

(1) 

(2) CP-violating ttj vertex 

iV2va j2 ft^ 5 /A 2 , (2) 

(3) CP-conserving jjH vertex 

-4vahi [ (kik 2 )g^ - k lv k^ }/A 2 , (3) 

(4) CP-violating jjH vertex 

Sva h2 k^e p<TllI/ /A 2 , (4) 

where v ~ 250 GeV, k and k\ i2 are incoming photon momenta, and a 7 i )7 2,Ai,fc2 are 
defined as 

a 7 i = sin6wR&{oc u w) + cos ^Re(a„ B ) , (5) 



a 7 2 = sin 8 w lm(a u w) + cos 6wIm(a' uB ) , (6) 
Qffti = sin 2 6 ) H/Re(a ¥ ,vi/) + cos 2 ^Re^s) — 2 sin#w cos 9 w^(<^wb), (7) 
0^2 = sin 2 9 w Re(a v y / ) + cos 2 ^ H /Re(a (/: ,^) - sin^ cos 8 w Re(a^ / B ) , (8) 

a, and a'- being the coefficients of Oi and O'j (i = uW, (pW, <pB, WB, (pW, ipB, 
WB and j = uB) respectively, and 9w the Weinberg angle. It will be helpful to 
note that the SM ffj coupling in our scheme is given by eQ/7 M , where e is the 
proton charge and Qf is /'s electric charge in e unit (e.g. Q u = 2/3). 

The top-quark decay vertex is also affected by some dim-6 operators. For 
the on-mass-shell W boson it will be sufficient to consider just the following tbW 
amplitude when nib is neglected: 



r wtb = u (Pb 



7*(1 -75) " V^/ftl + Tte) 



u{ Pt ), (9) 



where is given by 



f? = ~-^ M w[ ~^uw + ^ai Du , (10) 

with a Bu the coefficient of the operator Odu ■ On the other hand, the viW vertex 
is assumed to receive negligible contributions from physics beyond the SM. 

Finally, the initial-state polarizations are characterized by the initial electron 
and positron longitudinal polarizations P e and P g , the average helicities of the 
initial-laser photons P 7 and P 7 , and their maximum average linear polarizations Pt 
and Pi with the azimuthal angles ip and (p (defined in the same way as in 3\). The 
polarizations P T)t and $ have to satisfy 

0<P 2 + P t 2 <l, 0<P 2 + P^<1. (11) 

3. Optimal-observable analysis 

The calculation of the cross section is straightforward; to derive distributions of 
secondary fermions we have applied the Kawasaki-Shirafuji-Tsai technique JH] 
with FORM [2D] used for the necessary algebraic manipulations. We neglect con- 
tributions that are quadratic in non-standard interactions and treat the decaying 

"■"■Note that there is another potential source of contribution to f^, which may come from Oq u . 
However, this operator could be eliminated using equations of motion; therefore, it is neglected 
hereafter. We thank Ilya Ginzburg for pointing this to us. 



t and W as on-shell particles; therefore the angular-energy distribution of the b 
quark in the ee CM frame can be expressed as 

da 



fsu(E b , cos9 b ) + a 7 if 7 i(E b , cos 6^) + a y2 ME b , cos6 b ) 



dE b d cos 9b 

+ a hl f hl (E b , cos 9 b ) + a h2 fh2(E b ,cos9b) + a d f d (E b , cos 9 b ), (12) 

where fi(E b , cos 9b) are calculable functions; /sm denotes the standard-model contri- 
bution, / 7 i i7 2 describe, respectively, the anomalous CP-conserving and CP- violating 
tt7-vertices contributions, fhi,h2 those generated by the anomalous CP-conserving 
and CP-violating 77if-vertices, and fd that by the anomalous tbW-vertex with 

a d = Re(/*). 

Their analytical form is however too long to be presented in this letter. 

In order to apply the 00 technique, we first have to calculate the following 
matrix elements using the weighting functions fi(E b ,cos9 b ) defined in eq. (fT2"j) : 

Mij = J dE b dcos9 b fi(E b , cos 9 b )fj(E b , cos 9 b )/ f S M(E b , cos 9 b ), (13) 

and its inverse matrix Xy, where i, j = 1, ■ ■ • , 6 correspond to SM, 7I, 72, hi, h2 
and d respectively. Then, according to [21] , the expected statistical uncertainty for 
the measurements of a, is given by 



\Aoi\ = yJIoXa/Nt, (14) 

where 

I = J dE b d cos 9 h fsM (E b , cos 9 b ) 

and Nb is the total number of collected events. 

Inverting the matrix Ai, we have noticed that the numerical results for 
are often unstable: even a tiny variation of Aiij changes Xij significantly. This 
indicates that some of fi have similar shapes^ 2 and therefore their coefficients 
cannot be disentangled easily. Indeed, we already encountered a similar trouble 
in our latest analysis using final leptons ^H] . It is not surprising that we meet this 



" 2 Note that if two fi functions were proportional to each other, then the matrix would 
have a vanishing determinant, and therefore its inverse could not be determined. 



problem again here, since the main structure of the cross section is determined by 
that of 77 — > ti for both processes. 

The presence of such instability forces us to refrain from determining all the 
couplings at once through this process alone. Therefore, hereafter, we assume 
that some of a^s have been measured in other processes (e.g. in ee —>■ ti — ► £ ± X). 
Fortunately, however, we obtain some complementary information on coupling con- 
stants, which was not available in our previous analysis [T3], where only leptonic 
distributions were employed. 

Below we list all the elements of M (= M T ), which were computed for 

v /5^ = 500 GeV and A=l TeV. (15) 



(1) Linear polarization 

We chose the following values as typical linear polarizations: P e = P s = 1, P t = 
Pi = P 1 = P^ = l/y/2 and x{= V 9 — 0) = tt/4, where (p and (p are the azimuthal 
angles of P t and Pj. They are the same polarizations as those we used in [To] . 
1-1) m H = 100 GeV 



M n 


= 0.368 x 10 2 , 


M 12 


= 


.787 x 10 2 , 


M 13 


= -0.323 x 10 1 


M u 


= -0.145 x 10 2 , 


M 15 




0.153 x 10 1 , 


M w 


= 0, 


M22 


= 0.169 x 10 3 , 


M 23 




0.699 x 10 1 , 


M24 


= -0.299 x 10 2 


M25 


= -0.331 x 10 1 , 


M 26 


= 


.277 x 10 1 , 


M33 


= 0.352, 


M 34 


= 0.122 x 10 1 , 


M 35 


= 


.182, 


M3Q 


= -0.454, 


Mu 


= 0.681 x 10 1 , 


Ma 5 


= 


.583, 


M m 


= 0.271 x 10\ 


M 55 


= 0.987 x 10" 1 , 


M 56 




-0.281, 


M m 


= 0.866 x 10 1 . 


m H = 


= 300 GeV 












Mn 


= 0.368 x 10 2 , 


Mi 2 


= 


.787 x 10 2 , 


M Vi 


= -0.323 x 10 1 


M 1A 


= -0.359 x 10 2 , 


M 15 




■0.691 x 10 1 , 


M w 


= 0, 


M22 


= 0.169 x 10 3 , 


M 23 




0.699 x 10 1 , 


M 2 4 


= -0.742 x 10 2 


M25 


= -0.146 x 10 2 , 


M 2 6 


= 


.277 x 10 1 , 


M33 


= 0.352, 


M 3i 


= 0.298 x 10 1 , 


M 35 


= 


.681, 


M 3e 


= -0.454, 


M M 


= 0.421 x 10 2 , 


M& 


= 


.725 x 10 1 , 


Mag 


= 0.711 x 10 1 , 


M 55 


= 0.146 x 10\ 


M 56 


= 


.143, 


M m 


= 0.866 x 10 1 . 



(16) 



(17) 



1-3) m H = 500 GeV 



M u 


= 0.368 x 10 2 , 


M 12 


= 0.787 x 10 2 , 


M 13 


= -0.323 x 10 1 


M u 


= 0.170 x 10 2 , 


M 15 


= -0.101 x 10 2 , 


M w 


= 0, 


M22 


= 0.169 x 10 3 , 


M 23 


= -0.699 x 10 1 , 


M24 


= 0.352 x 10 2 , 


M25 


= -0.206 x 10 2 , 


M 2e 


= 0.277 x 10\ 


M 33 


= 0.352, 


M34 


= -0.148 x 10\ 


M 35 


= 0.809, 


M 36 


= -0.454, 


M44 


= 0.935 x 10\ 


Mm 


= -0.579 x 10\ 


M m 


= -0.283 x 10 1 


M 55 


= 0.369 x 10 1 , 


Mm 


= 0.253 x 10 1 , 


M m 


= 0.866 x 10 1 . 



;is) 



(2) Circular polarization 

We took the following values as circular-polarization parameters: P e = Pg = P 7 
= 1, which were also used in [T3j . 



2-1) m H = 


-- 100 GeV 












Mn - 


- 209 x 10 2 


KA 1 




n 454 x 1 2 

V . t: '.J Tfc s\ X \J . 


M iq - 
j vl 13 


- 


KA . . - 

J v 1 14 


- _n 690 x 10 1 


■> v 1 15 




—0 109 x 10~ 3 


J v 1 lb 


= 




- 988 x 10 2 


KA r,n 




n 


j vi 24 - 


- — n 144 x 10 2 

— \J . X XX s\ 1U 


KAnK = 


- —0 227 x 10 -3 


A/foc 
J vl 2b 




126 x 10 1 


Ml? = 


= 


- /vt 34 


- 


KA nt 




n 




- 


KAaa = 
J v l 44 


- 284 x 10 1 


AA ac 




457 x 10" 4 


•f vl m 


= 133 x 10 1 


A<55 = 


= 0.739 x 10~ 9 , 


Mm 




0.243 x 10~ 4 , 


M m -- 


= 0.393 x 10 1 . 


2-2) m H = 


= 300 GeV 












Mu 


= 0.209 x 10 2 , 


M12 




0.454 x 10 2 , 


M 13 = 


0. 


M u 


= -0.178 x 10 2 , 


M 15 




-0.177 x 10 1 , 


M 16 = 


0, 


M22 


= 0.988 x 10 2 , 


M 23 




0, 


M 2 4 = 


-0.373 x 10 2 , 


M 25 


= -0.368 x 10 1 , 


M 26 




0.126 x 10 1 , 


M33 = 


0, 


-M34 


= 0, 


M 35 




0, 


Mm = 


0, 


M44 


= 0.191 x 10 2 , 


Mm 




0.193 x 10 1 , 


Mm = 


0.360 x 10 1 , 


M 55 


= 0.198, 


M 56 




0.419, 


Mm = 


0.393 x 10 1 . 


2-3) m H = 


= 500 GeV 












Mu 


= 0.209 x 10 2 , 


M n 




0.454 x 10 2 , 


M 13 = 


0. 


M 14 


= 0.762 x 10 1 , 


M 15 




-0.502 x 10 1 , 


M w = 


0, 


M22 


= 0.988 x 10 2 , 


M 23 




0, 


M 2 i = 


0.159 x 10 2 , 


M 25 


= -0.105 x 10 2 , 


M 26 




0.126 x 10 1 , 


M33 = 


0, 


-M34 


= 0, 


M 35 




0, 


Mm = 


0, 


M44 


= 0.347 x 10 1 , 


Mm 




-0.233 x 10 1 , 


Mm = 


-0.138 x 10\ 


M 55 


= 0.158 x 10 1 , 


Mm 




0.103 x 10 1 , 


Mm = 


0.393 x 10 1 . 



(19) 



(20) 



(21) 



All the elements Mij above are given in units of fb. In these results, the third 
components of M for the circular polarization vanish This is common for 



analyses using the final lepton and the final 6-quark. Also, as in the leptonic case, 
-Mi6 = 0. This time, however, it is not because of the decoupling, which holds for 
the lepton production [22], but simply because the tbW vertex cannot contribute 
to the total cross section of 77 — > ti — > bX, since 



and Br(t — > bX) = 1, whatever anomalous terms are added to the tbW coupling 
as long as we assume that a top quark always decays through t — > bW. 

When estimating the statistical uncertainty in simultaneous measurements, e.g. 

of atji and a hi (assuming all other coefficients are known), we need only the com- 

[31 

ponents with indices 1, 2 and 4. Let us express the resultant uncertainties as Aa y { 
and Aa h {, where "3" shows that we used the input Adij, keeping three decimal 
places. In order to see how stable the results are, we also computed Aa I an d 
Aa h [ by rounding Aiij off to two decimal places. Then, if both of the deviations 

l^ a 7ifti — ^ a 7i hi\l ^°^\,hi are l ess t nan 10%, we accept the result as a stable 
solution. 

Although we did not find any stable solution in the three-parameter analysis, 
we did find some solutions in a two-parameter analysis; for those, the numerical 
results are presented below. According to the above criterion, the uncertainties for 
the following standard deviations Aoii are limited to 10%: 

1) Linear polarization 
• Independent of m# 



<x t ot(77 — tt -> bX) = Br(t -> bX)a tot (^ -> ti) 




(22) 



• m H = 100 GeV 




(23) 



• m H = 300 GeV 




(24) 



Aa hl = 5.4/ JN b , Aa d = 4.9/ 'JN b , (25) 



• m H = 500 GeV 



Aa l2 = 23/y/Nb, Aa h i = 5.0/ J N b , (26) 

Aa hl = 18/y/N b , Aa h2 = 22/^N b , (27) 

Aa hl = 8.0/y/N b , Aa d = 3.3//^, (28) 

where N b ~ 18400 for a luminosity of = eL eE = 500 ftT 1 with e being the 
relevant detection efficiency and L eS being the integrated luminosity" 3 

2) Circular polarization 

• m H = 100 GeV 

Aa hl = 14/ = 5.2/^, (29) 

• mir = 500 GeV 

A* M = 10/^, = 4.2/y^, (30) 

where N b ~ 10500 for Lf = 500 fb -1 . 

It is worth while to compare estimations of sensitivities obtained here for the bX 
final state, with those found in ^H] in the case of the final state. Unfortunately, 
here, we did not find any stable solution that would allow for a determination of 
<Xyi; the same was also observed for l^X. We therefore have to look for other 
suitable processes to determine this parameter. The precision of a 7 2 is not very 
good either, but it is still much better than in the case of the lepton analysis. On the 
other hand, we can see that analyzing the 6-quark process with linearly polarized 
beams enables us to estimate some Aa.i that were unstable in the lepton analysis, 
i.e. cases (J2H) and (|27j) . One of them, eq. (J2ZJ), is especially useful to probe the CP 
properties of heavy Higgs bosons through the determination of a hi and oth2- As for 
the determination of a^, the £ ± X final state seems to be more appropriate. These 
comparisons show that both final states (bX and £^X) provide complementary 
information and should therefore be included in a complete analysis. 



" 3 Hereafter we use the tree-level SM formula for computing Nb, therefore, below we have the 
same Nb for different m#. Also, for illustration, we assumed L eg — 500 fb _1 (adopting e = 1) as 
the standard reference point. However, one should not forget that tagging a 6-quark jet including 
its charge identification is harder than that of a lepton. 



The above results are for A = 1 TeV. When one takes the new-physics scale 
to be A' = XA, then all the above results (AaA are replaced with Z\ctj/A 2 , which 
means that the right-hand sides of eqs. p2j) - (JHUj) are multiplied by A 2 . 

4. Summary and Discussion 

We studied here beyond the SM effects in the process 77 — > tt —>■ bX for arbitrarily 
polarized photon beams, taking advantage of the fact that polarizations of the 
incoming-photon beams can be controlled. Non-standard interactions have been 
parameterized through dim-6 local and gauge-symmetric effective operators a la 
Buchmiiller and Wyler jni- Assuming that those new-physics effects are small, we 
have kept only terms linear in corrections to the SM tree-level vertices. 

We applied the optimal-observable technique to final 6-quark distributions, and 
estimated statistical significances of measuring each (allowed by the gauge invari- 
ance) non-standard parameter. Unfortunately, we had to conclude that it is never 
possible to determine all the independent non-standard parameters at once through 
77 — > tt — > bX alone. However, we still would be able to perform a useful analysis 
if we could utilize the complementary information collected in other independent 
processes. 

Comments on the background are here in order. The most serious background is 
iy-boson pair productions. Indeed, its total cross section could be 300 times larger 
than a tot (tt). Fortunately, however, a simulation study has shown that tt events 
can be selected with a signal-to-background ratio of 10 by imposing appropriate 
invariant-mass constraints on the final-particle momenta There, an efficient 
6-quark tagging is crucial, which is a basic assumption in the analysis presented 
here. 

Some non-standard couplings, which should be determined here, could also be 
studied in the standard e + e~ option of a linear collider. Therefore, it is worth while 
to compare the potential power of the two options. As far as the parameter a 7 i 
is concerned, the 77 collider does not allow for its determination, while it could 
be determined at e + e~. The second tt^y coupling a 7 2, which is proportional to 
the real part of the top-quark electric dipole moment,^ 4 can be measured here. It 
" 4 See [23 taking into account that the operators O u b, O q B and G q w are redundant. 



should be recalled that energy and polar-angle distributions of leptons and 6-quarks 
in e + e~ colliders are sensitive only to the imaginary part of the electric dipole 
moment," 5 while here the real part could be determined. For the measurement of 
couplings, e + e~ colliders are, of course, useless, while here, for the bX final 
state both at hi and could be measured. In the case of the decay form factor a rf 
measurement, the e + e~ option seems to be a little more advantageous, especially 
if e + e~ polarization can be tuned appropriately [25] . 

It should be emphasized here that the effective-operator strategy adopted in this 
article is valid only for A ^> v ~ 250 GeV, in contrast to the analysis of e + e~ — > 
tt — > t^X performed in [22j and [25] for example. Should the reaction 77 — > ti — > 
bX exhibit a deviation from the SM predictions that cannot be described properly 
within this framework, we would have an indication of low-energy beyond-the- 
SM physics, e.g. two-Higgs-doublet models with new scalar degrees of freedom of 
relatively low mass scale. 
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